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Abstract 

An incomplete particle identification distorts the observed event-by-event fluctuations of the hadron 
chemical composition in nucleus-nucleus collisions. A new experimental technique called the identity 
method was recently proposed. It eliminated the misidentification problem for one specific combination 
of the second moments in a system of two hadron species. In the present paper this method is extended 
to calculate all the second moments in a system with arbitrary number of hadron species. Special 
linear combinations of the second moments are introduced. These combinations are presented in terms 
of single-particle variables and can be found experimentally from the event-by-event averaging. The 
mathematical problem is then reduced to solving a system of linear equations. The effect of incomplete 
particle identification is fully eliminated from the final results. 
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I. INTRODUCTION 



A study of event-by-event (e-by-e) fluctuations in high-energy nucleus- nucleus (A+A) colli- 
sions opens new possibilities to investigate properties of strongly interacting matter (see, e.g., 
review [l| and references therein). Specific fluctuations can signal the onset of deconfinement 
when the collision energy becomes sufficiently high to create the quark-gluon plasma [2| at the 
initial stage of A+A collision. By measuring the fluctuations, one may also observe effects 
caused by dynamical instabilities when the expanding system goes through the 1-st order tran- 
sition line between the quark-eluon plasma and the hadron gas 13 1. Furthermore, the QCD 

Jlfl 

critical point may be signaled by a characteristic fluctuation pattern |4H6[. Fluctuations of the 
chemical (particle-type) composition of hadronic final states in A+A collisions are expected 
to be sensitive to the phase transition between hadronic and partonic matter. First data on 
the e-by-e chemical fluctuations from the CERN SPS |z-9| and BNL RHIC Q have been 
already published, and more systematic measurements are in progress. The e-by-e fluctuations 
of hadron multiplicities have been studied theoretically in statistical models (see, e.g., Ref. |ll| ) 
and in dynamical transport models (see, e.g., review 12( and references therein). 

Studies of the e-by-e chemical fluctuations assume particle number measurements for differ- 
ent hadron species (e.g., pions, kaons, and protons). The NA49 Collaboration jz-9| has used 
the measure o"d yn , which is defined as the difference between fluctuations observed in real and 
mixed events. The STAR Collaboration 101 has used, in addition to the aHvn measure, the 

nfi 



quantity z/ dyn (see, e.g., Ref. [13]). Moreover, it was suggested long ago [14|, ll5[ to quantify 



chemical fluctuations by the measure $ 



16| . Note that different fluctuation measures can 



be presented as specific combinations of the second moments of the multiplicity distribution. 
Some important features of different measures for the e-by-e fluctuations have been considered 



in Ref. 



12| 



A serious experimental problem of the e-by-e measurements of the chemical fluctuations is 
incomplete particle identification; that is the impossibility to determine uniquely the type of 
each detected particle. The effect of particle misidentifications distorts the measured fluctu- 
ations. For this reason the analysis of chemical fluctuations is usually performed in a small 
acceptance, where particle identification is relatively reliable. However, an important part of 
the information on e-by-e fluctuations in full phase space is then lost. Although it is usually im- 
possible to identify each detected particle, one can nevertheless determine with a high accuracy 



the average multiplicities (averaged over many events) for different hadron species. 



In Ref. 



18| a new experimental technique called the identity method was proposed. It solves 



the misidentification problem for one specific combination of the second moments in a system 
of two hadron species ('kaons' and 'pions'). In the present study we extend these results 
in two directions. First, we prove that not only the one specific combination of the second 
moments but all the second moments themselves can be uniquely reconstructed in spite of the 
effects of incomplete identification. Second, the identity method is extended to an arbitrary 
number k > 2 of hadron species. This is important for practical purposes since typically 
there is an incomplete identification for pions, kaons, and protons, which means k = 3. The 
identity method is introduced in Section [TTJ In Section II III the main results are presented. 
We also discuss several examples which illustrate some limiting cases of particle identification. 
Section IIVI presents the summary. 



II. IDENTITY METHOD 



The identity method was proposed in Ref. 



18] for the fluctuation measure This method is 



based on the fact that the analysis of chemical fluctuations can be performed within two different 
but fully equivalent formulations. The first formulation [l^] uses the identity variables; that is, 
the \l/ measure of chemical fluctuations is calculated using single-particle variables 
where the over-bar denotes averaging over the single-particle inclusive distribution. The event 
variable Z, which is a multi-particle analog of z, is defined as Z = J2iLi( x i ~ where the 
sum runs over the N particles in a given event. The measure ^ is defined as 

(N) ' 1 ' 

where the symbol (. . .) corresponds to the e-by-e averaging. One defines the single-particle 
variable Xj as the identity variable w±(i) which equals 1 if the ith particle is of the first type 
('kaon'), and W\(i) = if the ith particle is of the second type ('pion'). In a real measurement, 
it is unknown exactly whether a given particle is 'kaon' or 'pion'. As a consequence of this 
incomplete identification the variable Wi(i) is not exactly or 1, but becomes a distribution 
function with possible values in the whole [0, 1] interval. Nevertheless, despite the incomplete 
particle identification, one can directly use the definition (TTJ) to evaluate 
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In the second formulation, \I/ is calculated in terms of the moments of the multiplicity 



distribution. In the case of complete particle identification it was found 
1 



15| that 



(NY 



(Nt)(N 2 y + (NtfiN*) - 2(N 1 )(N 2 )(N 1 N 2 ) - (Ni) (N 2 ) - (Ni)(N 2 )' 



(2) 



where indices 1 and 2 correspond to different hadron species ('kaon' and 'pion'), and N = 
N\ + N 2 . Using the presentation (QQ), it was shown 18J that the measure \l/ can be factorized 
into a coefficient that represents the effect of misidentification, and the quantity (j2j), which 
corresponds to the value that ^ would have for complete identification. 



We follow Ref. [18] and assume that particle identification is achieved by measuring the parti- 
cle mass m. Since any measurement is of finite resolution, we deal with continuous distributions 
of observed masses denoted as Pj(m) and normalized as (j — 1, . . . , k > 2) 



dm Pj(m) = (Nj) . 



(3) 



Note that the functions Pj(m) are found for the different particle species using the values 
averaged over all particles from all collision events. The identity variables Wj(m) will be defined 

as 

k 



Wj (m) EE 



Pi' 



m) 



p{m) 



(4) 



i=l 



The complete identification (CI) of particles corresponds to distributions Pj{m) which do not 
overlap. In this case, Wj = for all particle species i ^ j and Wj = 1 for the jth species. When 
the distributions Pj{m) overlap, Wj(m) can take the value of any real number from [0, 1]. 
We introduce the quantities Wj, with j — 1, ■ ■ ■ ,k, and W p W q , with 1 < p < q < k, 

N(n) N(n) 
, W p W q EE ( W v( m i)) >< ( Yl 



N(n) 
i=l 



[mi 



(5) 



i=l 



i=l 



and define their event averages as 



(6) 



i -"cv 1 ev 

JVev n=l JVev n=l 

where A^ ev is the number of events, and N{n) = N\(n) + • ■ ■ + N^(n) is the total mul- 
tiplicity in the nth event. Each experimental event is characterized by a set of particle 
masses {mi,m 2 , . . . ,m^}, for which one can calculate the full sets of identity variables: 



{wj(mi),Wj(m2), ■ ■ ■ ,Wj(misr)}, with j = l,...,k. Thus, the quantities W 2 and W p W q are 
completely defined for each event, and their average values ([6]) can be found experimentally by 
straightforward e-by-e averaging. In the case of CI, one finds Wf = N 2 and W p W q = N p N q , 
thus, Eq. dSJ) yields 



(W 2 ) = (N 2 ) , (W p W q ) = (N p N q ) . 



(7) 



III. SECOND MOMENTS OF CHEMICAL FLUCTUATIONS 

The quantities (W 2 ) and (W q W p ) can be calculated as follows 

N 1= 0N2=0 N k =0 J J 

x J dm\P 2 (m 2 2 ) ...J dm k N2 P 2 (m 2 N2 ) x . . . x J dm k P k (m k ) ...J dm k Nk P k {m k Nk ) 

x \wj{m\) + ■■■ Wj{m l Ni ) + Wj{m\) + ■■■ + Wj(m 2 N2 ) + ... + w j (m'l) + • • • + ^(m^J 



(8) 



i=l 



Ki<Kfc 



OO OO LAJ /i n 

= EE-E • • • ' N *) / • • • / dm^P^ml 

N!=0N 2 =0 N k =0 J 

x J dm\P 2 (m 2 2 ) ...J dm k N2 P 2 {m 2 N2 ) x . . . x J dm k l P k {m k 1 ) ...J dm k Nk P } 



(ml) H w p (m} N ) + w p (ml) H h w p (m 2 N2 ) + . . . + w p (m k ) H h w p (m k N ) 



x w g (m}) H ^(mjvj + w q (m\) H h w q (m 2 N2 ) + . . . + w q (m k ) H h w q (m k N ) 



(9) 



i=l i=l l<i<l<k 

In Eqs. (jHJ) and (j^J), V(N%, . . . , iVfc) is the multiplicity distribution, Pi(m) = pi(m)/ (Ni) are the 
mass probability distributions of the ith species, and (s = 1,2) 

1 s 1 



dmw s Am) pi(m) , w ? 



dmw p (m)w q (m) pi(m) . (10) 



In the case of CI, when the distributions Pj(m) do not overlap, one finds that 



(11) 
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and Eqs. (jHJ) and ([9]) reduce then to Eq. (jZJ). The incomplete particle identification transforms 
the second moments (Nj) and (N p N q ) to the quantities (Wf) and (W p W q ), respectively. Each 
of the later quantities contains linear combinations of all the first and second moments, (iVj) 
and (Nf), as well as all the correlation terms (iVjiV/). 
Having introduced the notations 

k k 



( w f) - E w K - («i*) 2 ] = &i . « - E w [ 



i=l 



i=l 



one can transform Eqs. OH]) and dHJ) to the following form: 



E^ 2 K+ 2 E wviK 



j = l,2,...,fc 



i=l 

k 



Ki<Kk 



y](Nj) UpiU qi + (NiNl) [UpiUql + U p lU q ij = b pq , 1 <p < q < k . 



(13) 

— v , (14) 

i=l l<i<l<k 

The right-hand side of Eqs. ffl3|) and (fT4|) defined by Eq. f fl2|) are experimentally measurable 
quantities. The same is true for the coefficients (with s = 1 and 2) entering the left-hand 
side of Eqs. f lT3"|) and ( lTi|) . Therefore, Eqs. ( TT3|) and ( 1T3|) represent a system of + — l)/2 
linear equations for the second moments (N?) with j = 1, . . . , k and /c(fc — l)/2 correlators 
(N p Nq) with 1 < p < q < k. 

In order to solve Eqs. (JT3]) and flUj) we introduce the [fc + fc(/c - l)/2] x [fc + fc(fc - l)/2] 
matrix A 



( a\ 



A 



'12 



\ a l2 



(I 



(fe-l)fe 



.12 



(fc-l)fe\ 



.12 



a 



(jfc-l)* 



'12 



a 



(k-l)k 
12 



r/ 



12 

(fc-l)fc 



(15) 



a (k-l)k/ 
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where 



U.; 



a 



pq 



2ui P Ui q , 1 < p < q < k 



a lg = u pi u qi , 1 <p < q < k , i = l 
a 1 ™ = U p lU qm + U q lU pm , 1 < p < q < k , 1 <l <m < k 



1, . . . , k ; 



, k ; 



(16) 
(17) 
(18) 



The solution of Eqs. (TT3T) and (fl4|) can be presented by Cramer's formulas in terms of the 
determinants 

r2N det Aj , „ „ , det A r , 



det A 



*pg 
det A ' 



(19) 



where the matrices A~ and A, g are obtained by substituting in the matrix A the column 



J 3 3 3 

a \l ■ ■ ■ J Q fc> a i2> • • • 5 a (fc-l)fc 



and the column af 9 , . . . 



„P9 „P9 



> a fc 5 a 12> • • • > " (£_!)*:> 



respectively, for the col- 



umn bi, . . . , fefc, 6i2, • • • , b(k-i)k- Therefore, if detA ^ 0, the system of linear equations f JT3|) and 
pij) has a unique solution (fl9|) for all the second moments. In the case of CI (11 ip . one finds 
detA = 1, detA,- = 6j, and detA pg = 6 pr The solution (TIT?]) reduces then to Eq. ([7p. 
Introducing the [A; + — l)/2]-vectors 



( (Nf) \ 

(N 2 k ) 
(N 1 N 2 ) 



( h \ 



B = 



h 



(20) 



\(N k ^N k )J \b( k -i)kj 
one can write Eqs. ( 1X3]) and ( I14p in the matrix form AM = B. The solution ([19]) can be then 



rewritten as 



A/" 



A" 1 B 



(21) 



where A 1 is the inverse matrix of A. For two particle species, k — 2, this solution takes the 
form 



< (Nl) N 

(Ni) 

\(N!N 2 ) J 



( 



«21 



MliM 2 l 



«12 
M 22 



U\ 2 U 22 



2u n ui2 

2u 2 \U 22 
U U U 22 + M12M2I 



-1 



(22) 
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Then Eq. O yields 



/AT 2\ h l U 22 + h 2U 2 l2 ~ 2b 12 U 12 U 22 

\ iV l/ — , n2 ' l^" 3 / 

(_ttH«22 - U 12 U 2 i) 

, N 2s = Mi + ML - 2&12W11 ^ 

(Wll«22 - «12«2l) 

, , &12(miiM22 + M12M21) - biu 22 u 2 i - b 2 u n ui 2 

(iVi7V 2 ) - -2 ■ (25) 

{UuU 22 - M12M21J 

These results, inserted into Eq. fl2]), provide an alternative way to evaluate \1/ directly from the 
moments of the multiplicity distribution. 

In general, the particle- by-particle identification is difficult; that is, it is not known whether a 
given particle really corresponds to the j'th sort. On the other hand, the statistical identification 
in terms of the functions Pj{m) is usually reliable. Experimental measurements of the Pj{m) 
functions give the average numbers of each particle species. In most cases, a unique calculation 
of the second moments using Eq. ffT9l) is also possible. There is, however, an extreme situation 
when the only available experimental information consists of the average particle multiplicities. 
This leads to random identification (RI) which only defines, for each particle, the probabilities 
Pj of being of the jth sort. These probabilities are evidently equal to pj = (Nj)/(N), where 
iV = Yli=i N{- This situation is described by the mass distributions given by 

p.jim) = (AT.) f(m) , (26) 

where J dmf(m) = 1; that is, all functions Pj{m) have the same shape f(m) but different 
normalization (Nj). With these distributions one finds 

«Ji = W = Pj > j = l,...,k. (27) 

This leads to detA = 0, and Eqs. (TT3|) and (|T4|) do not define the second moments in a unique 
way. In fact, from Eq. (j5J) follows 

(Wf) = p) (N 2 ) , (W pq ) = p p p q (iV 2 ) , (28) 

that is, in the case of RI, the measured values (j28~|) include only the average multiplicities (Nj) 
and the second moment of the total multiplicity (iV 2 ). Equations (IT51) and (JHJ) for (N 2 ), and 
(NpN q ) reduce to a single relation 

k 

^(iV 2 ) + 2 (W = (iV 2 ) , (29) 

j=l l<P<9<fc 
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where the right-hand side of Eq. fj29|) is the only experimentally measured combination of the 
second moments. Therefore, RI gives only one restriction on k + k(k — l)/2 second moments 
and thus admits an infinite number of solutions for (AT 2 ) and (N p N g ). Any correctly normal- 
ized multiplicity distribution V{N i) . . . , N^), which reproduces experimental values of the first 
moments, would reproduce Eqs. ( TT3"]) and (THj) ; that is, in the case of RI the experimental data 
do not provide any non-trivial information on chemical fluctuations. 

It is also instructive to consider an illustrative example when particle species are divided 
into two groups: j = 1, . . . , k R with RI f )26|) and j — k R + 1, . . . , k with CI (JTTJ) . Equations f|T3|) 
and f )14p are then given by 

(K) = (wl) = ... = {w% R ) = (w 1 w 2 ) = ... = (w kR ^w kR ), n r = Y, N r, (so) 

(N R N q ) = (W x W q ) = ... = (W kR W q ) , k R + 1 < q < k ; (31) 



(iV 2 ) = (Wf) , j = k R + l,...,k; (32) 
(NpN q ) — (W p W q ) , j = k R +l,...,k; k R + l<p<q<k. (33) 

Quantities (Wf) and (W p W q ) can be measured experimentally by using their definitions ac- 
cording to Eq. For the particle species with RI, as follows from Eqs. ( 13"U|) and ( 13"T|) . the 
second moments of particle multiplicities include only their total multiplicity N R . Therefore, 
one knows all individual average multiplicities (Nj), but as far as chemical fluctuations are 
concerned, all particles in the RI group, 1 < j < k R , look undistinguishable. On the other 
hand, this fact does not prevent calculations of the second moments ( )32|) and (|33|) in the CI 
group, k R + 1 < j < k. 

In the formulation considered in this paper, the functions Pj(m) are defined as quantities 
averaged over all particles from all collision events. One can consider the set of events with fixed 
total multiplicity N. All formulae of this paper straightforwardly apply in this case too. The 
only modifications are: 1) event averaging (. . .) over all events is changed to averaging (. . .)n 
over events with fixed A^; 2) the functions Pj(m) should be replaced by Pj(m; N) calculated for 
fixed N. This procedure may open some new possibilities in the studies of chemical fluctuations. 
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IV. SUMMARY 



An incomplete particle identification prevents a straightforward measurement of the second 



moments (N?) and (N p N q ) 
method proposed in Ref. 



the multiplicity distribution. In this paper we extend the identity 
. We introduce the quantities (Wj) and (W p W q ). Each of these 
quantities is a specific linear combination of all the first and second moments (Ni) and (Nf), as 
well as the correlation terms (NiNi). The quantities (W?) and (W p W q ) are presented in terms 
of e-by-e averages of functions depending on the single-particle identity variables according to 
Eq. (jSJ), and can thus be measured experimentally. Mathematically, the problem of finding all 
the second moments (N?) and (N p N q ) is then reduced to solving the system of k + k(k — l)/2 
linear equations (I13p and fll4p . All coefficients entering the left-hand side of these equations 
are given in terms of experimentally measurable density functions Pj{m). The right-hand side 
in Eqs. (TT3]) and (TT4]) is defined by Eq. (j!2p which also includes experimentally measurable 
quantities. In most cases the determinant of matrix (1151) is not equal to zero and, therefore, 
all second moments of particle number distributions can be uniquely reconstructed by Eq. (TTUj) 
from event-by-event measurements despite the effects of incomplete identification. This is 
valid for an arbitrary number k > 2 of different hadron species. The matrix A -1 in Eq. (|2ip 
represents the correction of the measured values f|T2l) . Such a correction eliminates the effect 
of misidentification. This provides the values of all the second moments (N?) and (N p N q ) in 
a model-independent way, as they would be obtained in an experiment in which each particle 
is uniquely identified. However, all measured quantities entering Eqs. (IT5|) and (THE]) contain 
experimental errors. Therefore, the practical applicability of the identity method procedure 
constructed in this paper requires further studies. 
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